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1 Introduction 



One of the characteristic features of quantum field theory is the non- conservation of the 
number of particles. The most familiar cases are the collision processes, where genuinely 
quantum mechanical interactions lead to creation (or annihilation) of particles. Here, all 
the interacting fields are quantum fields. It has been known for long that also classical 
external fields can change the number of particles. This is the well-known Schwinger 
mechanism for particle production in an external field IQ. The Schwinger formula has 
been derived in a number of ways and has found a wide range of applications, from 
Hawking radiation to heavy ion collisions [Q. 

Schwinger's result has some shortcomings. The external field is homogeneous through- 
out space; the timescale is large (infinite) and it applies only at zero temperature. Further- 
more the probability of finding the system in the groundstate after some time has elapsed 
is not related in a simple manner to observable quantities, like the average number of 
particles created because the probability distribution is non-trivial. 

The first problem has to some extent been discussed elsewhere and will not be 
considered here. 

The problem of large times comes about because of the nature of the model, which 
consists of free particles minimally coupled to an external background field. In this case 
the backreaction, i.e. the field generated by created particles is neglected. Since the num- 
ber of particles, and hence the induced field, increases with time there is a maximal time 
for which the background field approximation applies. If one insists on that approxima- 
tion, which we do, it is important to keep the explicit time dependence throughout the 
calculations. The backreaction has been studied extensively in 

The third problem of extending the pair creation phenomenon to non-zero temper- 
atures has been studied in thermal field theory, where the effective action is a natural 
object of interest. However, there seems to be confusion about the influence of a thermal 
background on the pair creation P, 0, §]. In the real time formalism one does not find a 
thermal correction whereas using imaginary time there appears to be one. It is also not 
clear if any of the formalisms actually are applicable to this kind of problems. In this 
context it is again necessary to keep the time dependence and care must be taken with 
the infinite time limit. We find that at fixed temperature, the thermal contribution to 
pair production will vanish in the infinite time limit. At any finite time, however, there 
is a significant thermal effect. In the background field approximation, one can therefore 
question the use of effective lagrangians (the infinite time limit being understood) at least 
concerning the imaginary part responsible for pair production. 

So much for the physics, we now turn to the formalism. We will be working in the 
functional Schrodinger representation (see e.g. [TD|, [TT|) which is convenient for many 



reasons. Because of the close analogy with quantum mechanics, there are no conceptual 
problems in this picture. We do not suffer from doubts about the formalism and/or appli- 
cability, as in thermal field theory. Furthermore, the treatment of time and temperature 
is straightforward. 

In order to make the calculations explicit we make extensive use of the so called "field 
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basis" concept, i.e. one views the wave functionals as overlaps of a general state and a 
"field basis state". This is known to be possible for bosons. We extend the formalism to 
fermions. One then has to be careful with the Grassmann nature of the basis. Also extra 
care must be taken when working in a reducible representation, which is the case at hand. 

With this setup, the strategy is simple. We solve the time evolution for the thermal 
density matrix of the system under consideration, i.e. we find p satisfying 

if)=[H,p\ , p(0) = equilibrium ensemble. (1) 

Then we simply calculate the thermal average of the number operator A^, that is 

{N) = iv [N p{t)]. (2) 

Closely related to this quantity is the number density in momentum space (transverse 
and longitudinal momenta). We also find the probability distribution of the n-particle 
states in the ensemble at arbitrary times, specifically we calculate the probability of 
finding the system in the groundstate at finite times. The paper consists of two main 
sections. In the first section we treat bosons and set up the machinery. Fermions are 
dealt with in the second section. After a careful discussion of the fermionic field basis, 
all the calculations previously done for bosons are repeated for fermions. The paper ends 
with a few concluding remarks. 

A word of caution; throughout this paper there are numerous seemingly infinite quan- 
tities. This is not case however, since we always work with a UV-cutoff. Only in the 
expressions that are not only convergent but also quite insensitive to the cutoff, do we 
remove this cutoff. 

2 Bosons 

The hamiltonian for a charged scalar field in an external electromagnetic field is 

H = j £xt({^)tt\-k) + j rf^a;rf^#^(x)cj2(x,y)0(y), (3) 

where 

^2(x, y) = (-(V. - zeA(x))2 + m2)5(x - y) (4) 

and we use the convention that e = — |e| < 0. In a standard manner ||^, ^ we represent 
the field operators and the states: 

^^(j) , TT ^ , 0^ ^ (/)* , Tft ^ (5) 

t d(p I 0(p 

, ^vl/(0,0*)*. (6) 

By * we mean complex conjugation. The inner product is defined by functional integration 

= (7) 
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where = VcpVcj)* . As a reference we give a formula for gaussian functional integrals: 



V'^cp exp 



d^x<fy (f)*n(f) + j {a*(f) + (j)*b) 

= det{n)-^ exp J d^xd^y a*VL-^h 



where the determinant is defined through det VL = exp Tr In Q with a functional logarithm 
and Tt Q = J d^x r2(x, x). 

The state functional can be viewed as an overlap between a general state and a "field 
basis state" , thus 

vI/(0,0*) = (0,0*|M/). (9) 



The field basis has the following properties 
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For operators we have 



tr O 

(0,0*|d(0,0t^7r,7rt)|vl/) 



5(0-V^)5(</>*-^*), 



V'<h(A<J)*\0\ 



1, 



(10) 
(11) 
(12) 



(13) 
(14) 

(15) 



2.1 The initial ensemble 

The system under consideration is described by a density matrix p{t). For a good intro- 
duction to density matrices see |T^. We intend to evolve the system in time, using the 
hamiltonian (^. The system evolves from an initial ensemble at thermal equilibrium at 
t = to a final ensemble ai t = tj when the external field is switched off and the number 
of particles in the system is counted. 

Initially the system is in thermal equilibrium and is described by an unnormalized 



stationary density matrix p„ = e ^ satisfying 

— dppu = HqPu 



(16) 



where Hq is the hamiltonian (|^) with A = and P is the inverse temperature, P = 
Note that such an ensemble cannot be created by pure QED since it contains all possible 
charges. Furthermore we have 



limpn = 1 , Pi = Pu- 



(17) 
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In the functional representation p„ is represented by a functional matrix 0202) — 

(0i,0t|p„|02,0;), thus 

S 5 

- dppuiMl, 0202) = -^o(01, 01, T:^' TTi)Pu{Mh </'202), (18) 

501 d0J 

and the two additional conditions read: 

limp„(0i0t,020;) = (5(0i-02)(5(0t-0;) (19) 

P«(010I,0202) = P^(0202,010l)- (20) 

In order to solve for pu, we make a gaussian ansatz with a covariance ft^j 



-/d3^d3^0*(x)Q°.(x,y)0,(y) 



i,j = l,2. (21) 



p„(010*,0202) = -?V"«exp 

One obtains 

- df) In = Tr 1]?^ (22) 

-dgQ'l, = O?i'-c^o' (23) 

-a^0?2 = (24) 

~t^/3^21 = ^21^11 (25) 

-dpQl, = (26) 
where loq — u; having set A = 0. We have also used the notation 

(m)(x,y) = J dhQ{^,z)Q{z,y). (27) 

Due to translational invariancc the covarianccs will satisfy Q" (x, y) = ^ij(x — y). We 
can therefore simplify the equations by performing a fourier transform according to 

0(x) = / 0(p)e*P '' , 7r(x) = / 7r(p)e-^P-^ , 0(x - y) = / Q{p)e'^<--y^ (28) 

where Jp — J ci^p/(27r)^. The equations then read 

-a^ln7V„ = V [ 0?i(p) (29) 
Jp 

-a^Q?i(p) = Q?i(p)'-a;o(p)' (30) 

-df^n^ip) = Q?i(p)Q?2(p) (31) 

-df^nUp) = ^^^i(p)^^?i(p) (32) 

-d^nlM = (33) 

where a;o(p) = \/p^~+ln^ and V — (27r)^5^(0) is the large but finite spatial volume. The 
hermiticity conditions on yield 

^?1(P)* = ^22(P) , ^^?2(P)* = ^?2(P) , ^^2l(P)* = ^^2l(P) (34) 
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The solution to (|29|)-(p^) leading to a density matrix that satisfies the delta function 
condition (|19]) is 



fi?i(p) 

^?2(P) 



iV„. 



^22 (p) — t^o coth /Jcuo 
sinh /JtUo 



^21 (P) 



exp 



V / In 



sinh Pluq 



The precise value of A^^ can be found using lim/3^0 / T^'^'PiPu 
from (plQl) . The partition function Z = tr and the normalized density matrix 
Po = Pw/tr p„ are now readily found as 



(35) 
(36) 

(37) 

1 which follows 



Z = exp 

Po(010I,0202 

exp 



-V I 21n(2sinh(^)) 

V / ln(2c(Jotanh 
Jp 



(3J 



- / ((0^01 + 0202)t^ocoth/3c^o 



X 



sinh (3ijiJQ 



(0t02 + 02</'l)) 



(39) 



2.2 The time dependent ensemble 

We solve the Liouville equation for the density matrix p(t) in an external field A(t) with 
the initial condition specified in the previous section, i.e. 



idtp{t) = [H{t)rp{t)] 

P(0) = Po. 

The external field is constant in space but varies in time as 



(40) 
(41) 



,t < 
E(t) = { Ee^ ,0<t<tf A{t) 
,tf<t 



,t<0 
-Ete^ ,0<t<tf 

-EtfBs ,tf <t 



(42) 



where i? > so that eE < 0. In the functional representation we find the equation for 

P(010t,0202): 

s s s s 

idtp{(pi4>l, Ml) = H{(j)i, 0*, — , T--)p(0i0t, 0202) - H{(j)2, (j)*2, J—, t-t)p(0i01' 0202) 

001 d(pi 0(p2 0(p2 

(43) 

Again we make a gaussian ansatz, 



p(0i0t,020;,t) = Ar(t)exp 



- j d?xd^y 0*(x)%(x,y,t)0,(y) 



(44) 
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and perform a fourier transform as in the previos case (with the specified external field 
the hamiltonian is still translationally invariant). Suppressing the time and momentum 
dependence, the equations obtained are 

idtlnN = V f{nu-n22) (45) 
Jp 

idtn^ = - - ^]l2^^2l (46) 

tdtni2 = ^]l2(^^ll-^^22) (47) 

tdtn2l = l]2l(f^ll -^^22) (48) 
^a^^22 = -r^22' + - ^^12^^21 (49) 



where uj{p) = y (p — eA)^ + m^. The density matrix p is hermitian and the hermiticity 
conditions read 

fill* = r222 ) fil2* = fil2 ) ^21* = ^21 (50) 

By inspection one finds that the solution to three of the equations are 

fil2 = a(fiii+fi22) (51) 
^21 = 6(^11+^22) (52) 
^22 = fill* (53) 

where a and b are constants of motion determined by the initial conditions. From (|35| ) and 
([36| ) we find a = b = — 1/2 cosh /3u;o. Eliminating fii2,fi2i and fi22 we find the remaining 
two equations 



idt In 



N = 1/ /(fi - fi*) (54) 
Jp 

idtfl = fi^-cu^ \ (fi + fi*)2 (55) 

4 cosh i3uJo 

where we have set fin = fi. Hence the problem has been reduced to finding fi. Now write 

Vt = —i-^ y = r exp {iO coth (3uJo) . (56) 

y 

Inserting these expressions into (|55|) and rewriting the initial condition, (l35|), we find 

,-,2 



„2/ 



f - -f + cuV = , r(0) = 1 r(0) = 



r'e = ujo , ^(0) = 6(0) = ooo (57) 
Here we have introduced dimensionless quantities 



t = t^^J\eE\, p = p/y\eE\, m = m/^\eE\, (3 = (3^\eE\, 
A = (pi)2 + (p2)2 ^ ^2^ ^ _ U + ^pi _ t^2^ V = V\eE\'^/^ 



(5J 
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and the overdots in ( |5^ ) denote d/dt. Finally introducing the complex function 
z = rexp{i9) the set of equations (E7^ can be collected into the form 



z + u z 
z{0) 

m 




1 



(59) 



One can now find a closed expression for Q in terms of special functions, since the general 
solution of (BSf) is a linear combination of two parabolic cylinder functions [O 



D 



2'^ 2 



±(l+0(t-p3)). 



(60) 



The solution satisfying the initial conditions is easily written down but since the expression 
is rather long we don't give it here. Instead we end this section with the expression for 
the time dependent normalized density matrix, 



p{t) = exp 



exp 



V J InUn + n*) tanh Pcoq tanh 



X 



1 



2 cosh Puq 



where Q in terms of r 



1 2; I reads 
Q{t) = 



/ UJQ COth/^LUp _ T 



(61) 



(62) 



2.3 Expectation value of the number operator 

Once the density matrix has been found one can start asking questions about the sys- 
tem. We will be concerned with particle creation and therefore concentrate on thermal 
expectation values of the number operator. 



{N){t) = tr (iVp(t)) 

= J V^{D^26{ct>, - 02)5(01 - 0DiV(0i, 0Dp(0i0t, 020;; t). (63) 

More specifically we will calculate the number of particles at t = tf just after the field has 
been turned off. We then find ourselves in the gauge A = —Etje^. In order to make use 
of operators expressed in terms of fields in the A = gauge we make a transformation to 
that gauge. The fields transform according to 

A(a;) ^ A{x) - VX{x) 0(a;) (j){x)e-''^^''\ (64) 

Since the physical states are gauge invariant the covariance transforms as 

Q{x, y) ^ fi5(x, y) = Q{x, y)e~^^(m-m) _ (65) 
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For the specific transformation tliat clianges A = —Etfe^, into A = 0, tlie effect on tlie 
covariance in tlie density matrix ( 16^ simply amounts to a shift in the momentum 

(66) 

(67) 
(68) 
(69) 



P ^ Pp = P + eA{tf) = {p\p\p'^ - eEtf) 
Hence r,uj and c^o transform as follows 



c.o(p) - ujl{p) = Jk + {p^ + ifY 



Now find the number operator N . The gauge invariant stationary groundstate, ex- 
pressed in terms of fields in the A = gauge, is given by 



^^{(pAl = det(cJo)"'/'exp 
Associated with this groundstate are the annihilation operators 

^ (a;o(p)0(p)+?7rt(p)) 

(^o(p)0^(p) +^7r(p)). 



a[p) 



Hp) 



2uo{p) 
1 

2cc;o(p) 



(70) 

(71) 
(72) 



The states created by a^(p) have charge e and physical momentum p, while states created 
by &^(p) have charge (— e) and physical momentum (— p), as seen from the momentum 
and charge operators expressed in terms of a's and 6's, 



p(at(p)a(p)-6t(p)6(p)) 



Q = ie I d'x (0^7r^ -(j)n) = e I (aT(p)a(p) - ¥ip)b{p)) 



The expression for the number operator in the A = gauge is 



N 



;at(p)a(p)+6t(p)6(p)) 
' ^ vrt(p)vr(p)+c.o(p)0^(p)0(p)-V-). 



p UJo[P) 



(73) 
(74) 

(75) 
(76) 



By (^3D and the gauge transformed versions of (|6TD and (^2]), we find after some manip- 
ulations the result 



V 
2Jp 



-,9 



rt + -^r., ^) coth 



. ^0 2 I 

^ r,9 ^9 ' ,-, r,9 '9 ' ,-, '9 



r,9 



UJo^Q 



(77) 
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Figure 1: The momentum distribution, n(— p)/2, of the positively charged bosons at zero 
temperature. All quantities are expressed in units of the electric field. 



To find the number of particles actually created during the process (-/V)", we subtract 
the number of particles at t = 0: 



{N)\t=o = V 



coth ^—^ — 1 



Hence, we obtain 

{Ny'{tf) ^ V I n(p,i/)coth 
Jp 



\^ roo roo 

dp^ dAn{A,p\tf)coth 



(78) 



(79) 



where we have defined the zero temperature number density of the created particles 
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Figure 2: The average number of created bosons for rh = 1. The dashed line shows the 
asymptotic behaviour at zero temperature. All quantities are expressed in units of the 
electric field. 

n{p,tf) as 

Equations (|79|) and (^) constitute the exact expressions for the expectation value of the 
numbers of particles created by the external field at finite temperature, when a time tf 
has elapsed. Since, as we have argued, the positively charged particles created by &^(p) 
have a momentum (— p), their momentum distribution at zero temperature will be given 
by n(— p)/2 = (&^(— p)&(— p))/3=oo- Figure |I| shows this distribution. Note that for small 
times, there are very few particles with zero momentum. The result of integrating (|79D 
numerically for a mass m = 1 is shown in figure |^. We see that increasing the temperature 
increases the number of created particles but for large enough tf the slope is independent 
of temperature. 

We have so far only considered a simple time dependence of the electromagnetic field, 
but it is straightforward to generalize to an arbitrary time dependence. Introduce the two 
dimensionless functions f{t), F(t) and let 

r ,t < 

E(t) = I Ef{t)e^ ,0<t<tf Ait) = 
[ ,tf <t 

The functions f{t) and F{t) must satisfy 

F{i) = f{i) , F(0) = (82) 



,t < 

-EF{t)/.J\^\ eg ,Q<t<tf (81) 

-EF{tf)/^f\^\ 63 ,tf<t 



10 



As a matter of fact, (^) and ( pOD still hold, with the modifications that Vg satisfies the 
gauge transformed version of (|57D, where now 

= ^Ja + (p3 - F{t) + F{if)Y (83) 

and similarily for ujq and uJq. Because of the new time dependence, however, cannot 
be expressed in terms of special functions, but has to be found numerically. This is in 
practice not a serious drawback, since the integral still must be evaluated numerically. 



2.4 Asymptotic expansions 

We have obtained the exact expression for the number of particles created for arbitrary 
times tf. It is, however, possible to obtain some more convenient expressions when t/ is 
large. To that end we investigate the behaviour of the number density n. The number 
density may be written 

„(A, p^ i,) = ^fs)\'+\i,m _ 1 (84) 

'''''I' u,„(zl{t,)z,{t,)-z,{t,)z-^(t,)) 
where Zg{t) by (|59D satisfies 

■ig + {t-if-p'fZg = Q. (85) 

Since the number density n is not sensitive to a rescaling of Zg, the relevant initial 
condition for Zg following from ( ^9] ) is 

%(0) 



%(0) 



iul (86) 



Let us now be more specific and consider a large final time tj » 1. For such t/, consider 
A and satisfying 



-tj - VA « -VA , p-^«-l , » VA. (87) 

For A and p^ restricted by (pTl) the initial condition (|86D reads 

i{if+p^ 



%(0) 

A solution to (^) satisfying ([88| ) is 

zg{i) = D_,^_, (-(1 + 0(t - tf - p')) , (89) 
which may be confirmed using the first of the following asymptotic expansions of parabolic 



cylinder functions when |x| >> >> ||T3 



D^{xe'-^) ~ e-'^ixe'-^y , x>0 (90) 



, ZTT . ix'^ . iTT , , , '\/ '^iTT ix^ ^ V , . 1 , 

DJxe-) ~ e-—ixe-Y-—, -e" e"*''^ xe" " , x < 0. 91 

F -z/ 
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It is now a straightforward calculation to find the asymptotic behaviour of Zgiif) in which 
case (|91|) applies. By (^), and using u)q ~ \p^\ and |r(|A + \)\~'^ = \ cosh one finds 
the number density 

r^(A,p^^/) ~ 2e~^^. (92) 

What happens if p lies outside the region in (|H7|)? It is clear that for \p^\ sufficiently large 
compared to if, we can neglect the time dependence in the equation for Zg. Hence there 
will be no created particles i.e. n (A,p^^/) = 0. In the region where |p| is neither large, 
nor restricted by (|87|), the contribution to {N)" will be small compared to the dominant 
contribution. Therefore we can effectively set n to zero in this region. We conclude that 
for tf » 1 and if » rh, 



2e-^^, -tf <p^ <0, 

uyiv,p ,if) = \ 

and thus 



0, otherwise. 



V fv poo ill}'' 

m) = ^ "f i 2e- coth (94) 
Two limits involving (3 may now be considered, namely 1/(3 « if and 1/(3 » if, 

{Nr{if) = -^.ife--"^ = -VEHfe-^ , ^ « if (95) 



(Arr(t7) = —\n{2if)e-^ ={-)y'^\n[2tf^\eE\je-l^, ^ » if m 



The result ( |95D is identical to the zero temperature result obtained in [Q, eq. (5.1.61). 



2.5 Probabilities 

The probability -P(^E') of finding the system in a normalized state |\E') is 

Pi^) = (97) 

We are interested in the probability P{n,m;tf ) of finding the system in a state with n 
particles and m anti particles at time t = tf. For specific momenta these states are 

n m 

|p,...p., k,...k^) = n n a\p.)b\k,)\0)- (98) 

The desired probability is a sum of probabilities for such states with n and m held fixed, 
i.e. 

P{n,m;tf) =V''-^'^ / (pi...p„, ki...k„|p(t/)|pi...p„, ki...k„) (99) 

One finds that the probabilities, P{n,m;tf), are generated by a function f{x,y]tf) as 



Qn+m 

P(,v,n;«,) = /(i,M/)-;^^;^/(..y;*/) 



(100) 

x'=i;=0 
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It immediately follows that J2^m=oP{^^''^'^tf) — 1- The generating function f{x,y;tf) 
reads 

fix, y;tj)= exp [-Tr In {(1 - xg{tf)){l - yg{tj)) - xyh{tf)}] (101) 

where the functions g{tf) and h{tf) are expressed in terms of the number density n{tf) 
defined in ([80|). Suppressing the momentum dependence, we have 



h{tf) 



1 



(n(t/) + 1 + coth pu^) sinh Plj^ 



[n 



(102) 
(103) 



We note that when j3 = oo only h{tf) survives. Since h{tf) multiplies xy we conclude 
that only uncharged states are present at zero temperature. This is what we expect since 
all states have been created out of the uncharged groundstate. On the other hand when 
t/ = we have that h{0) = while g{0) = exp(— /^tuo), leading to ordinary Boltzmann 
probabilities for the initial ensemble. 

The probability P(0, 0; tf) is readily obtained 



exp 



- coth 

2 ^ 2 



— Tr In 



n 



{tf) + l + coth/3( 



(104) 



At zero temperature this reduces to 
P(0,0;t/) = exp 
which for large times, tf » 1, by (^ 
P(0,0;t/) = exp 



-Tr ln(l + ^) 



becomes 
Vt 



13 = oo 



(105) 



2(27r 



■ttA\ 



(106) 



This is the result obtained by Schwinger in 1951 |T|. We end this section by noting that 
expectation values of some operators may be obtained from the generating function. For 
the number operator we have 



which of course agrees with (ffTl) . 



[n + 'm)P{n,m;tf) 

n,m=0 

f{l,l;tj)-'d.J{x,x;tj)l^^ 



Tr [{n{tf) + l)coth 



(107) 
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3 Fermions 

We will set up the formalism by studying an illustrative example, the two-dimensional 
fermionic harmonic oscillator. The annihilation and creation operators, a,b and d\b\ 
satisfy 

{a,a^} = {6,6^} = l. (108) 

All other anticommutators vanish. The Hilbert space is spanned by the four orthogonal 
states |0,0),|1,0),|0, 1) and |1,1), 

a|0,0) = 6|0,0) = 0, a^|0,0) = |1,0), 6^|0, 0) = |0, 1), a^6^|0, 0) = |1, 1). (109) 

The states are normalized, (0,0|0,0) = (1,0|1,0) = (0, 1|0, 1) = (1,1|1,1) = 1- We 
represent this system by functions ^(77*, 77) of a complex Grassmann number 77 and its 
complex conjugate 77* in the following manner. 



|o,o)^i, 11,0)^77*, 10,1)^77, 11,1)^77*77, 

a-«r^drf*, a^^-^T)*, b d^, b^ rj, 



(110) 



so that overlaps with Grassmann states arc given by 

(77*7710, 0) = 1, (77*7711,0) =77*, (77*7710,1) = 77, (77*7711,1) =77*77 (111) 

and 

1 

{V*,VW*,V') = E {V*,v\n,m){n,m\r]'*,r]') (112) 

n,m=0 

= l + ri*ri' -ri'*ri + ri*riri'*ri' ^e'^*'''-''"'^. (113) 

The Grassmann state may now be expanded as 
1 

|,7*,7)= J2 |n,7n)(7i,7n|77*77) = |0,0) + 11,0)77 +10,1)77* + 11,1)77*77. (114) 

n,m=0 

The partition of unity can thus be expressed as 

i = |o,o)(o,o| + |i,o)(i,o| + |o,i)(o,i| + |i,i)(i,i| 

= j (fr]'(fr]\r]*r]){r]*r]\r]'*r]'){r]'*r]'\, (115) 
where (Prj = dri*dri. Hence an arbitrary overlap is given by 

(*l|*2) = / dVA(*l|^*r7)(77 * 77 |77'*77')(77'Vl*2) 

= ydVAe"*"'""'*" *t(7^*,7^)*2(V*,V)- (116) 
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Let O be an operator, 

1 

0= Yl C>nin2,mim2\ni,n2){mi,m2\. (117) 

ni ,n2 ,mi ,m2=0 

The trace of this operator is then, 

tr a = Ooo,oo + (^10,10 + Ooifii + OuM- (118) 
This can be expressed as, 

ti O = f d^r]'d'^T]{T]*r]\d\T]'*r]'){T]*r]\r]'*T]'). (119) 



The "naive" trace, 

J d'r]'d'r]{r]'*r]'\d\v*7^){r]*7^W*r]') = 0oo,oo - ^lo.io - Ooi.oi + (120) 

is in fact the supertrace Str((9) = tT{0{—l)^). Let us now turn to field theory. The 
Hamiltonian for Dirac fermions in an external electromagnetic field A is, 

H = J dVy^['4jl{^),My)]ho^p{^,y), (121) 

where 

h{x, y) = (-Z7% ■ (V. - 2eA(x)) + m7°)5(x - y). (122) 
The anticommutation relations read 

{M^)Jl{y)} = 5^p5{^-y). (123) 

All other anticommutators vanish. As was first done in [P], we let the fermionic field 
operators act on wavefunctionals "^{r]*, rf) of a complex Grassmann field t] and its complex 
conjugate rj* according to 

4(x) » -L(,„„ + _i_,, (124) 
^:(x) " ^(,:(x) + ^). (125) 



The representation defined by ( |124| ) and ( |125| ) of the field operator algebra is reducible. 
This can be seen explicitly by introducing the operators, 

1 A 
6'„(x) ^ ^(?7^(x 



^l(x) ^ ;^(-<(x) + ^). (126) 
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We have 

{i(x),^;(y)} = 5,^5(x-y), (127) 

and 9 as well as 9'^ anticommute with both ip and ip'^ making the representation reducible. 
Contrary to what is said in we may now view wavefunctionals as overlaps with Grass- 
mann field states I?]*//), 



^{r]*,r]) = {T]*r]\^). (128) 
The partition of unity is given by functional Grassmann integration, 

1 = J D^r]'D^r]\r]*r]){r]*r]\r]'*r]'){r]'*r]'\, (129) 

having set D^t] = Dri*Dri and hence the inner product is 

(^i|^2) = / D^D^T] {V*VW*V') ^l{v\v)^2iv'*,v')- (130) 

For an operator O we have 

{rj*rj\0\ri'*ri') = 0{jfrj,rj'*rj') 
(r/*r/|a(^,V^t)|vi,) = a(r/*,^,r/,|-)(r^*r/|vl/), 

trrf(C)) = j D^7]'D^r]{r]*7]\d\r]'*7]'){7]*r]\r]'*7]'). (131) 

Since the representation is reducible, all states will have an unphysical degeneracy with 
the degree of degeneracy given by d. Hence the physical trace is given by 

tr(a) = ^tr,(a). (132) 

The degeneracy d is given by the number of states generated by 9"^ and 9. For symmetry 
reasons this is equal to the number of states generated by ip"^ and ip. The total number 
of states in the functional space is tr(i(l). Hence, 

d^ = trrf(i) = J D^ri'D^ri{ri*ri\ri'*riy = det(4/) = exp(4y ln4), (133) 
where /a/3(x,y) = SapS^yi - y). Thus 

d = det(2J) = exp{AV / in 2). (134) 

Jp 

Let us end this subsection by giving a standard formula for Gaussian functional integrals, 
D^r]exp{ J d^xd'^y r]l{x)n^^{x,y)r]fi{y) + J d^x <(x)xa(x) + C(x)r/a(x)) = 

det{-n) exp(- J d'xdS CMn-'p{^, y)xpiy)) (135) 



16 



3.1 Initial ensemble 



We will now do all the calculations for fermions that were previously done for bosons. Let 
Ho denote the Hamiltonian with A = 0. By (|16D and (|l^) we find, letting 

Pu{vlVl,V2V2) = {V*l,Vl\pu\V2,V2), 



d 

limp„(77*r/i,r/2*?/2) 



6 6 



(136) 
(137) 
(138) 



{vhVi\V2^V2), 

Pu{V*lVuV2V2) = P*u{V*2V2,V*lVl)- 

Make a Gaussian ansatz (suppressing spinor indices), 

Puiv*iVuV*2V2) = iV.exp( J d'xd^y r/*(x)fiO.(x, y)r7,(y)) t,j = 1,2. (139) 
Using (|136|) we obtain 



5/3lniV. 



-Tr(/i°fi?i), 



(/ - r^?i)/i°(/ + fi? 



(/ - fi°i)/i°fi° 



125 

■ 



—9/3^22 



1 



2^21^ ^12' 

where is h having set A = 0. Furthermore, by ( p.37| ) 

1, 



lim 
lim 
lim n° 



hm 1^22 = 0, 
- lim = 



and by (PSI) 



^22) (^12)^ ~ ^12' (^21)^ ~ ^ 



21) 



ipx 



where (f2"'')Q,/3(x, y) = fi^^(y,x). Fourier transforming, just like for bosons, 

fi(x,y) = / ^](p)e^P■(--y), r?(x) = / r/(p)e 
we find the solution satisfying ( |14(]|) - ([148 

= ^22(P)^ 



7°^^-P + "^^anh(^), 



fi?2(p) = fill(p) + l, fi^i(p)=^^?i(p)-l. 



exp 



/ In cosh ■ 



(140) 
(141) 
(142) 
(143) 
(144) 

(145) 
(146) 
(147) 

(148) 
(149) 

(150) 
(151) 
(152) 
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4V I In 2 cosh 
p 2 



Hence one may obtain the partition function 

Z = tipu = exp 
and the normahzed density matrix, po = Pu/trp„, 

Po{vlVi,V2V2) =exp{-4V / ln2)x 

Jp 



(153) 



exp 



VP 



(154) 



having suppressed the momentum dependence. 
3.2 The time dependent ensemble 

Let us now find p{t) with the initial condition p(0) = po- By (^0]), one has 



6 6 6 6 

idtp{v*iVuV2V2) = H{ril, — , r/i, ^)p(r/i*r/i, r/2^2) - (^(^2, j^, V2, ^)/^(^2^2, ^1*^))*- 

(155) 



We make an ansatz of the same form as the initial condition, 

p{v*iVi,V2V2,t) =exp{-AV [ ln2)x 

Jp 



exp 



((^r + V2) tanh(^) n{t) {rii + rj^) + rilrj^ - rilrji] 
p z 



(156) 



in momentum space. Since p is hermitian, fi(p, t) is a hermitian matrix. Hence, with the 
same external field as for bosons, using (|15(]|) we obtain 



idtn=[h,^l n(p,0) = -7° 



(7 ■ p + m) 
^o(p) 



Introducing the spinors Xi ^-nd X2, 



where, 

we solve ( |157| ). Noting that. 



1 + = XiXi + X2xL 
idtXj = hxj, i = l,2, 

^(p) = 7°(7 ■ (P - eA) + m) 



(157) 

(158) 
(159) 



and hence /i(p)^ = i^(p)^ = (p ~ eA)^ + m^, we differentiate ( |159D with respect to t and 
obtain, 

^ + {co' + 2eE7V)x, = 0. (160) 



To have ( |16CI| ) diagonal we choose the chiral representation of the gamma matrices, 



7 



-1 
-1 



, 7 



(T 

-a- 





-a^ 



(161) 



(162) 



in which, 

<T ■ (p — eA) —m 

—m —<7 ■ (p — eA) 

It remains to choose initial conditions for xi and X2 such that the initial condition for 
is satisfied. This is most easily done by diagonalizing fi(t = 0). Let 



^{^) = UDU\ x,(0) = f/xf, (163) 
for some unitary matrix U chosen in such a way that D = diag(l, 1, —1, —1). Explicitly: 



U 



m 





r+ 



m 



(164) 



where r_|_ = y2t<;o('^o + P^) and r_ = ^2ujq{u}o — p^). Thus, 

X?{x?y + X2{X2y = l + D = diag(2, 2, 0, 0), 



and we choose. 



X? 



V2 



1 






X2 



^-72 





1 






(165) 



(166) 



leading to. 



Xi(0) 



X2(0) 



UJo{uJo + p^) 



Uo{uo + p^) 



m 


p-^ +UJo 

p^ + ip^ 
—p^ + ip^ 

P^ + tuo 

m 



(167) 



168) 



Now hP is proportional to fi(0), = — ti;o^(0). This implies that both Xi(0) and X2(0) 
are eigenspinors of with the same eigenvalue {—uo). By ( |159| ) we then have. 



(169) 



19 



Hence we obtain, using ( p,6CI| ), 



1 



Xi 



n3^ 



X2 



myi 


{p^ + uo)y2 

{p^ + ip'^)yi 

{-p^ + ip'^)yi 
{p^ + uJo)y2 

rfiyi 



(170) 



(171) 



where, 



y\ + {u'^ - i)yi = 0, ?/i(0) = 1, ?/i(0) = iuo, 
y2 + (cu^ + i)y2 = 0, 1/2(0) = 1, 1/2(0) = iu^o, 



(172) 
(173) 



having introduced dimensionless quantities as for bosons, yi can be written as a hnear 
combination of the (hnearly dependent) parabohc cyhnder functions. 



D 



(-(l+z)r), D.^_,{{1+^)T), D.^{il-^)r), D.^{-{1-z)t), (174) 



whereas ?/2 can be written in terms of. 



D 



Note that, 



;-(l + z)r), D.^{{1 + z)t), D.^_,{il-f)T), D 



lA-l 



~(l-t)T). 



y2 



jp^ - i)yi - iyi 



(175) 
(176) 



3.3 Expectation value of the number operator 

The number operator, expressed in terms of fields in the A = gauge, is simply, 

Jp 1 ^Qy?) 



(177) 



As for bosons, we transform the covariance to the gauge A = by shifting the momentum, 

fi3(p,t/) = l](p + eA(tj),t^). (178) 
The expectation value of the number operator then follows. 



(iV)(t^) = tr(iV/3(t^)) = \l[^- tanh(MP) 
Subtracting the number of particles present at t = 0, 



)tr(^](p,o)fi3(p,^/)) . 



(iV)|,=o = 2r /(l-tanh(^)). 



(179) 



;i80) 
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Figure 3: The momentum distribution, n{—p)/2, of the positrons at zero temperature. All 
quantities are expressed in units of the electric field. 



the number of particles created is, 

(Nritf) = ^ /^tanh(M)(4-tr(Q'^Q^)). (181) 

Writing, 

= X?(XS)^ + X^(X^)^ - 1, Xiixiy + xlixlV - 1, (182) 

one finds, 

ir{QPn^)= Y: |(x°,xf)r-4, (183) 
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Figure 4: The average number of created fermions for m = 1. The dashed line shows 
the asymptotic behaviour at zero temperature. All quantities are expressed in units of the 
electric field. 



where (x, A) = Xa^a ^^"^ having used that the norm x) is conserved in time. Using 
( |17CI| ) and ( |171| ) we obtain, 



— |A2/? + (cuo + P^)(cUo^+P^)2/fi' 



Finally, writing yi = re*^ and again using the conservation of the spinor norm, we find 



;i84) 



{Nr{tf) = V I n{p,tf)tanh{^), 



having defined the zero temperature number density n{p,tf) as 

n(p,t/) = 2 

where 



+ uj'^r 



0, r(0) = 1, r(0) = 0, 
r', /(0) = cuo, 

4' ^(o) = o. 



(185) 

(186) 

(187) 
(188) 
(189) 



The subscript (superscript) g still indicates gauge transformed quantities i.e. p is shifted, 
pI = P^ + if. Figure ^ shows n(— p,t/)/2, the average momentum distribution of the 
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created positrons at zero temperature. Comparing with the corresponding resuh for 
bosons, figure there is an interesting difference for small times. Lots of fermions are 
created with zero momentum while hardly any bosons are. The result of numerically 
evaluating the integral ( |185| ) with a mass rh = 1 for some different temperatures is shown 
in figure |. Increasing the temperature clearly decreases the particle production but 
again, as for bosons, for large enough times the slope is independent of temperature. 
This suppression of particle production can be understood in a simple model at fixed 
spin and momentum. We then have four states: |0), |e~), |e"^) and \e~e~^). Due to charge 
conservation the only allowed particle production- and annihilation processes are: |0) — 
\e~e~^) , \e~e~^) — >■ |0). At infinite temperature the probabilities for |0) and \e~e~^) 
are equal in the initial ensemble. Hence creation and annihilation processes are equally 
probable leading to no creation on the average. 

If we have an electric field with an arbitrary time dependence as for bosons, (0), we 
have to change uj accordingly, gauge transformations become Pg = + F{if) and we also 



have to change (|188|) into / = f{t)r 



3.4 Asymptotic expansion 

Let us now examine what happens to the number density n{A,p^,if) when t/ is large, 
tf » 1. As for bosons, consider A and p^ satisfying. 



^A « tf, 
The initial condition for is, 

n«(p,0) 



A-tf « p-^ « 



p^ « -1. 



-7 



o(7-(p + eA(t/))+m) 



a;o(p + eA(tj)) 
Using (|90|) one finds ^^(p, t) for p restricted by (|190|) as. 



I — — ttA 

V2e— 



I — — ttA 

V2e— 





D_.^(-(l + z)(t-t;-p3)) 
i±i(pi + ^f)D_,^^_,{-(\ + ^){t-tf - f)) 

r l±i(-pi + ^p^)D_,^^_,{-{\ + - tf - f)) 
D.^{-{l + z)it-tf-f)) 



m 



D.^_,{-{l+z)(t-tf-f)) 



(190) 



(191) 



(192) 



(193) 



Note that the Xj here are not normalized exactly in the same way as in ( |1671 ) and ( |168D , 
the normalization differs by some irrelevant phases. Again for p restricted by (|190| ), we 
obtain, using (pl|). 



n(A,p^^/) = 2 - hiin'^ns (tf)) = 4e-"^. 



(194) 



23 



Thus, using a similar reasoning as for bosons, we conclude that for >> 1 and tf >> m, 

MA,P^t;)= , .Z_L (195) 



Comparing to the analogous expression for bosons, (|93D , we note that the only difference 
is a factor of 2 arising from the spin degrees of freedom of the fermions. Accordingly we 
obtain 

{Nntf) = dp' j^^ dA 4tanh(^)e-^ tf » 1, if » m. (196) 



Considering the limits ~ «tf and ~ >> tf of ( |196| ), we find 



— vrm^ TWIT'S J 



(27r)3 ' 7r2 ' ' /J 

(iV)-(t,) = ^t>--^ = (^)3/W|i?l^tJe-T^, (19J 



3.5 Probabilities 



We will now as for bosons give the expressions for the probability P{n, m; tf) of finding n 
electrons and m positrons in the ensemble at time t = tf. The calculations are however 
a good deal more involved for fermions so we will illustrate what happens by calculat- 
ing P{0,0;tf) explicitly. A particular choice of groundstate is, (remember that the 
representation is reducible), 

MV\V) = det(-^)exp(^r7*fi%). (199) 

We find in momentum space 

= det(^) det(fi° + n^) X 



det ((fi^ - i)(fi° + - 1)-' - {nf^ - + 1)(^]° + nf^y' 

det 



■^^^ nnf^ - i)(i]° + n^){nP + 1)(^]° + n^) ^i-th^^ 



A{th^ - 1) 4 



(200) 



where ^2^^ = th Q3{tf), th = tanh(^-^) and det^p denotes the determinant evaluated only 
with respect to the spinor indices. Now express and ^^{tf) in terms of spinors, see 
182|) , and use 



detsp{A) = ^[{trA)^ - 6{tiAY{tiA^) + 8trA tiA' + 3{tiAY - 6trA^], (201) 
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which holds for a 4 x 4-matrix A. Using a computer to handle the algebra we obtain 



P{0,0;tf) = {^i^,\p{tf)\^,) = det{-{l + th{\{xlxi)\'-2)+thy) 

= exp(2V / ln[i(l + th{2 - n{tf)) + th^)]), 
Jp 4 



(202) 



expressed in terms of the zero temperature number density n{tf) defined in ( |186|) . At 
zero temperature this reduces to 



P(0, 0; tf) = exp 



2Tr ln(l 



which for large times, t/ >> 1 , by ( |195|) becomes 



P(0,0;t/) = exp 



Vtj 



roo 

/ c/Aln(l-e-"^ 



(203) 



(204) 



Again this result agrees with the result of Schwinger In general we express the 

probability P{n,m;tf) in terms of a generating function f{x,y;tf) as for bosons, ( |100| ), 

fix, y- tf) = exp(2Tr ln[(l + x g{tf)){l + y g{tf)) + xy h{tf)]), 

where gitf) and hitf) are 

2 



(205) 



hitf) 



(2 - n(t j) + 2 coth f3Ql) sinh (3u)l ' 

n{tf){A-n(tf)) 
{2~n{tf) + 2coth/5i 



(206) 



Finally calculating the expectation value of the number operator 

oo 

{N){tf) = ^(n + m)P(n,m) = /(l,l)-ia,/(a;,x)U= 

n,m 

= Tr (2 + tanh(M)(^t;)_2)), 



(207) 



which agrees with ( |179| ). 



4 Final remarks 

We will briefly comment on the infinite time limit of P{0,0]tf) since this quantity is 
usually discussed. We will do this for fermions. For bosons, a similar argument may be 



made. The expression ( P02| ) for P(0,0,t/) may be rewritten when tf » 1. Using (|195|) 
we find for an arbitrary temperature 



P(0,0;t/) = exp 



C{P) + 



2V 



2 1-2 



oo pO 

dA 



, o , 1 + th(2 - 4e-"^) + th"^ 

ITTWp 



(208) 
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where C{(3) is a time independent constant and th = tanh(-^). The only time dependence 
is in u!Q{tf). From this expression it is clear that w, in the literature misleadingly known 
as the pair production rate, is 

„ = Ito !Ef(MM = ^r,A _ e--) (209) 

tf-*co Vtf TT Jrffi 

for any temperature. Therefore in the strict infinite time limit w is independent of tem- 
perature, even though for intermediate times it certainly does depend on temperature. 
However, the quantity that is directly related to the observable particle production is not 
w but the average of the number operator, (N)^"^, which we have thoroughly discussed. 
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